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Abstract. One of the objectives of theories describing quantum dynamical geometry is to 
compute expectation values of geometrical observables. The results of such computations can 
be affected by whether or not matter is taken into account. It is thus important to understand 
to what extent and to what effect matter can affect dynamical geometries. Using a simple 
model, it is shown that matter can effectively mold a geometry into an isotropic configuration. 
Implications for "atomistic" models of quantum geometry are briefly discussed. 



1. Introduction 

In the context of quantum gravity (see e.g. [1] for an overview of different approaches), much 
attention is devoted to studies of pure-gravity systems in which all matter degrees of freedom 
and their interactions are switched off. Whereas the use of this assumption is understandable 
given that most proposals for quantum theory of gravity are difficult to study even without the 
inclusion of additional fields, it is known that the presence of matter can alter the behavior 
of a gravitational system [2]. The purpose of this paper is to illustrate the effect of matter 
on geometry in a simple setting. In particular, the purpose is to show that matter can induce 
isotropy in a model that does not assume it a-priori - an effect not easily seen from standard 
treatments of matter in dynamical geometries [2] . 

Physical theories are often formulated using the Feynman path integral and many candidates 
for a theory of quantum gravity [1] are also defined within this framework. These candidates 
propose to study path integrals (partition functions) formally written as 

Z* = J V[g] j V[4>]e- E ^ (1) 

where integrals denote summations over geometries g and matter degrees of freedom cp, and 
E(g, (j)) is a weighting function that depends on both the geometry and the matter content (this 
function can be interpreted as the action or the energy of the system, according to the setting). 
The weighting function is often assumed to be of the form E(g, cp) = E g (g) + E ( p(g, <f>) so that 
the contribution from the matter fields in E ( p(g, (ft) is separated from the pure- geometry part 
E g (g). Furthermore, it is often assumed that if the matter contribution E^g, (j)) is set to zero, 
the resultant, simpler, path integral 

^=o = j V[g] e~ E ^ (2) 

has many of the same properties as the original one (1). That is, it is assumed that expectation 
values {O)^ and (0)^=0 of a geometric observable O computed in the two models have similar 



properties. It is this assumption that motivates the study of "pure-geometry" systems even 
when the real universe is known to contain matter fields of various types. 

The purpose of this paper is to demonstrate that geometric observables can in fact be 
different depending on whether or not matter is included in the calculation. The strategy 
to demonstrate this will be to compute expectation values explicitly in the setting of a two 
dimensional rectangular box in which the sides lengths L x and L y are variable but the total area 
L x L y = A is fixed. This system was already introduced in [3] but it is explored here more fully. 
An important property of this system is that it contains a single configuration where L x = L y 
and the box is a square but it contains a multitude of other configurations in which L x ^ L y ; 
the geometric observable R = L x /L y differentiates the square configuration from the non-square 
ones. The main result will be that the expectation value for this ratio in a system without 
matter must be different from one, but that it can be one when matter is present. 

The rest of the paper is organized as follows. Sec. 2 describes the flexible box system and 
defines its sum over geometries. The expectation value of the observable in the absence of matter 
is also computed. Sec. 3 discusses the effect of adding particles. It is shown through numerical 
studies that the expectation value of the observable R = L x /L y in the presence of particles can 
be close to unity. This indicates that matter can select a homogenous and isotropic geometry 
from the ensemble. In Sec. 4 this discussion is extended to the case of fields and it shown in 
what circumstances the Casimir energy of fields has a similar effect. A discussion of applications 
to current research programs on "atomistic" models 1 of spacetime are presented in Sec. 5. 



2. Toy Dynamical Geometries 

2.1. Regularized Sum Over Geometries 

Formal expressions like (1) and (2) do not specify the space of geometries to be integrated over or 
how exactly to carry out the integrals. To overcome this, several approaches (e.g. [3-5]) propose 
to study systems like (1) in a discretized setting. In those approaches, one replaces the integral 
over geometries by a finite sum over a set C of geometry configurations. The same strategy 
is adopted here. Since the purpose of this paper is to study effects of matter on dynamical 
geometries and not to propose a sensible way or regularizing the path integral for quantum 
gravity, the set C is chosen to make this specific application possible. 

Consider configurations spaces Ca corresponding to sets describing rectangular boxes having 
sides of length L x and L y and area A. To make each of these sets Ca finite, the possible lengths 
L x and L y must be restricted to a discrete or semi-discrete spectrum characterized by at least 
one length scale, £, designating a smallest allowed length. If the allowed lengths are assumed to 
be equally spaced when they are small, a concrete length spectrum that can be used is 

T ( W if L < A 

Le \R XL>\. (3) 

The minimal value and next larger values for L are £ and integer multiples of £, respectively, up 
to a scale A. Above this scale, L can take any real number value. It will be convenient to choose 
A = \J~A~. With this convention, the sets Ca can be written out explicitly as 

^-{(^),(«,^),....(V2,V3),...,(4,/)}. (4) 

Note that these sets are countable and well defined for all real A > £ 2 . 



1 The term "atomistic" model is due to H-T. Elze. 



2.2. Partition Function and Geometric Observables 

After defining the configuration space of a dynamical box with area A, one must specify how 
the configurations should be weighted in the ensemble defining the partition function (1). In 
the case where they are all equally weighted, the partition function of the system is 



x/e 

£<» 

P =i 



(5) 



The first term shown in the sum represents values of L x in the discrete part of the length 
spectrum, i.e. the first half of the configurations listed in (4). The second term is the analogous 
sum over discrete values of L y 2 . In the end, Z turns out to be equal to the size of the set Ca] 
the thermodynamic properties of the empty box system are determined entirely by entropy. 

Expectation values of observables describing the box geometry can now be computed. For 
concreteness, consider the ratio R = L x /L y . The expectation value for R should be close to 
unity if the box is likely to be square, and very different from unity if the box is likely to be 
long and narrow. One finds 
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(6) 



The last expression denotes the evaluation of (R) in the limit A — > oo. Since the result is 
divergent in this limit, (also when I — > 0), the calculation suggests that the box should be 
expected to be long and narrow in that limit. Curiously, note that the expectation value for the 
inverse ratio, (L y /L x ), is given by exactly the same diverging expression; this is a consequence 
of the L x <-> L y symmetry of the system and implies that these ratios provide measures for the 
degree to which the two sides of the box are different and cannot be used to infer which length, 
L x or L y , is larger. 



Another useful expectation value is 
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the last expression again stands for the A — ► oo limit. This quantity can be used to compute 
the fluctuations (Ai?) 2 = (R 2 ) — (R) 2 of the ratio R. In the limit A — > oo, the contribution to 
(AR) 2 from the first term, (R 2 ), is much larger than that from the second, (R) 2 , and thus the 
fluctuations diverge as A 3 / 2 . 



3. Particles in Dynamical Geometries 

3.1. Particles as Standing Waves 

Suppose that the flexible box is populated with quantum mechanical particles of mass m. The 
dynamics of each particle is described by the Schrodinger equation in two dimensions with a 
potential that is zero inside the box and infinity everywhere else. After setting H = c = 1, the 
energy spectrum for each particle is 

2 The given prescription for evaluating Z carries a risk of over-counting whenever the two configurations 
corresponding to p = \/£ are in fact equivalent. This problem can be averted by considering A such that 
X/t is not an integer. In any case, the error involved becomes negligible when the area is large. 



where k x and k y are the momenta in the x and y directions, respectively. The allowed values 
for these momenta (i = x,y) are 

h = (9) 
where a depends on the type of boundary conditions, 

1 Dirichlct 



Periodic, 

and rii = a, 1, 2, . . . are quantum numbers whose range also depends on the boundary conditions. 
For a given box configuration (L x ,L y ), the partition function for a single particle is 

z(L x ,L y ) = £ exp^ ^-^J E H 2^-^ (") 

where /? = 1/&T is the inverse temperature. For a system of iV identical Maxwell-Boltzmann 
particles, their combined contribution is 

z N (L x ,L y ) = (z(L x ,L y )) N /N\ (12) 

The total partition for the particles as well as the box can then be written as 

Z = 2j2^[pt,- e j (13) 

where the factor of two in front has the same origin as in (5). 

It is important to emphasize here that z and zn depend explicitly on the lengths of the sides 
of the box and the parameter a. In particular, it can be shown that the largest term in Z is due 
to configurations with L x = L y when a = 1, or to configurations with L x 3> L y and L y 3> L x 
when a = 0. These features form the basis for understanding the results described next. 



3.2. Evaluation of Observables 

The expectation value for R can now be computed, analogously to (6), using 

<*> = ^K<pO' + (14) 

This can be evaluated numerically given certain values for the parameters. To make for a realistic 
situation, the mass of the particles and the temperature of the system should be much smaller 
than the inverse minimal length, to, kT <C Thus /3(2 — a) 2 n 2 /8m 3> l/l 2 . After fixing this 
ratio in terms of £, which can then be set to unity without loss of generality, the remaining free 
variables are A and N. In the following, the number of particles iV is expressed in the form 
N = Nq + cA a with No and a positive dimensionless constants, and c a positive constant with 
dimension depending on the value of a. 

Numerical results for (R) and (Ai?) 2 with £ = 1, 0(2 — a) 2 '"it 2 /8m = 10 are shown in Fig. 1. 
Curves corresponding to a fixed number of particles (N = 1, a = 0, 1) show that the expectation 
value (R) and its fluctuations rise as the area of the boxes is increased. The asymptotic 
slopes match the behavior expected from the empty box calculations in (6) and (7). Curves 
corresponding to systems wherein the number of particles scales with area (N = 0.01^4, a = 0, 1) 




exhibit more interesting behavior. When Dirichlet boundary conditions are imposed (a = 1), 
the asymptotic behavior for large A is (R) — > 1 and ((AR) 2 ) — ► 0. Thus, in this case, the 
ensemble of particles mold the dynamical geometry into a stable isotropic configuration. The 
situation is different for periodic boundary conditions (a = 0); there both the expectation value 
for R and its fluctations diverge faster than in the empty box situations. 

The most interesting result of these plots is the emergence of isotropy, (R) — * 1, in the 
case of Dirichlet boundary conditions. It can be understood qualitatively as follows. The 
isotropic configuration contributes the most out of all the terms composing the partition function 
z(L x ,L y ) of a single particle. This advantage is then raised to a power cA a . Since the number 
of non- isotropic configurations grows at most polynomially with A (it grows as \f~A) , all of those 
configurations contribute negligibly to the overall partition function in the large area limit. The 
difference in growth rates between the polynomial and the exponential components also explains 
why the approach to perfect isotropy occurs rather quickly. 

4. Casimir Energy in Dynamical Geometries 

4-1- Regularization of Casimir Energy 

Another way of implementing matter is using fields instead of particles. The expectation value 
of the Hamiltonian operator of a field is in general given by a sum of two terms: one term 
proportional to the number of particles in some state (zero for the vacuum), and one other term 
that is divergent even in the absence of physical particles. This latter term is called the Casimir 
energy and it is usually attributed to quantum vacuum fluctuations [7]. It is (setting H = c = 1) 

i 

The factor rj is positive for bosonic fields and negative for fermionic fields. It is different from 
one if a field (such as a vector field) can have multiple polarizations. The summation is over all 
field mode energies. 

In a flat 2 + 1 dimensional background, relativistic field modes have energies uj given by 
to 2 = m 2 + k 2 + k 2 in terms of momenta k XjV and mass m. The momenta are again constrained 
by boundary conditions given in (9) and (10). Thus the Casimir energy is 

1 0000/ 2 2 \ I/ 2 

E = \wi2-a) £ £ ^ + || + ^ , (16) 

n x =a-n y =cr \ x V J 



where m CT = m/(2 — a)ir. This quantity is ultra-violet divergent; much work on the Casimir 
energy revolves around extracting its physically relevant components [7]. One technique to do 
this is to define a modified version of the summation (16) which includes a cutoff function Dg 
multiplying each summand as follows: 



oo oo 



£ ^1^L ^ D \t'L y ■ (17) 

n x =an y =a K y/ 

The cutoff function Dg must be chosen to satisfy certain criteria [8]. First, it must reduce to 
unity when £ — > so that Eq^ — > Eq in this limit. Second, for £ ^ 0, the function should fall off 
quickly when its arguments become large as to make Eq g finite. Third, it should not have any 
singularities or branch cuts in the complex plane. Lastly, it is also convenient to choose Dg so 
that it is symmetric in its arguments, n x /L x <-> n y /L y . 
To evaluate (17), one applies the Abel-Plana formula 

J>(n)=/ F(t)dt+-F(0) + i —-^(FW-Fi-it)) (18) 

n=0 J ° J ° 

twice in order to exchange the two summations over n x and n y into integrals. The resultant Eq^ 
can be written (see [7-9] for more details) in the form 

E 0/ = X - r)7r (2 - a) (E t>A + E e>B + Ear). (19) 
The first two terms in the parenthesis are 

Ee,A = (L x L y ) / du x / du y J ml + u 2 x + ul D e (u x ,u y ) oc 

Jo Jo v * (2Q) 

E €)B = - (1 - 2a) (L x + L„) ^ du ^/mJT^ D e (u,0) oc (1 - 2a) y . 

These terms are divergent when £ — > and thus parametrize the ultra-violet behavior of £?o- 
When £ ^ and m <C 1, their leading behavior is as shown. The precise coefficients are 
dependent on the details of the cutoff function. Interestingly, the terms are proportional to the 
area and perimeter of the two-dimensional box, and the latter has different sign for the two 
considered boundary conditions. 

The last term in (19), Ear, consists of the remaining contributions obtained from the 
expansions using the Abel-Plana formula. For the present discussion, the important features of 
Ear are that it is finite in the £ — > limit and that its magnitude decreases with area [6, 7, 9]. 

The Casimir energy is an active subject of research [6, 7] because it is a purely quantum 
effect that leads to a macroscopic effect in the form of the Casimir effect/force. In the standard 
description of the Casimir effect, it is the finite term Ear that is physically interesting. A clear 
explanation of why the divergent terms can be discarded when describing situations relevant for 
experiments testing the Casimir force is given in works on cavities with a piston [8]. In short, the 
reason is that the divergent contributions from inside and outside a cavity cancel from formulae 
describing measurable quantities (the force on the piston). 

In the context considered in this paper, the field is defined only within the boundaries of 
the box, and thus a cancellation of divergent terms from inside and outside cannot occur. The 
analysis of the effect of Casimir energy on the statistical properties of the dynamical box must 
therefore either include the regularized divergent terms, or discard them in an ad-hoc manner. 



4-2. Casimir Energy and Expected Geometry 

The partition function of the flexible box with field is 



VA/t 

Z = 2£^0M (21) 
P =i 

where z<f,{(3,p) stands for the partition function of a field theory at inverse temperature (3 in a 
box configuration labelled by p. should have contributions from the Casimir energy and 

also from the real particle present in a thermal state. However, in the low temperature regime 
where the average thermal energy is smaller than the mass of the particles, it is reasonable to 
approximate by the contribution from the Casimir energy alone. Thus z ( / > (P,p) ~ e _ ^ E °^ 
with Eq given by (19) for each box configuration labelled by p. 

Consider first the option of including the divergent but regularized terms of Eq in the 
calculations. Since the area of the boxes is kept fixed, the term Eg a contributing to Eq is 
the same for all the box configurations and so can be factored out and essentially ignored. The 
other terms depend on the geometry of the box and thus must all be taken into consideration. 
Being divergent in the £ — > limit, the term EgB proportional to the perimeter of the box is 
dominant over the term Ear- The latter term can therefore be ignored. The partition function 
thus becomes 

x/e , (AW 
Z = 2 £ exp (-I3 W {2 - a )(l-2a)p(pe + -)). (22) 

Expectation values can be computed numerically similarly as done in the previous section. The 
results are that the square configuration become dominant if the parameters a and r/ are chosen 
as either a = 1 and r/ < 0, or as a = and r] > 0. That is, the square configuration is preferred 
if either the field is fermionic with Dirichlet boundary conditions or if the field is bosonic with 
periodic boundary conditions. 

Next, consider the option of discarding the regularized divergent terms of Eq and using only 
the finite contribution Ear- Since Ear is decreasing with A, it's effect would be insignificant 
in the partition function for the dynamical box. The statistical properties of the dynamical 
box should thus be expected to match those described in Sec. 2. In other words, the finite 
component of the Casimir energy is not sufficient for shaping the dynamical geometry into a 
particular preferred configuration. 



5. Discussion 

Any successful theory of quantum dynamical geometry must eventually include coupling to the 
matter content observed in the universe. A-priori, average properties of dynamical geometries 
should be expected to be different when computed in theories with and without matter. 
Consequently, it is important to understand in detail the way in which matter may affect 
geometric observables. This issue was addressed in this paper using a toy model in which 
dynamical geometry was described by a two-dimensional "flexible" rectangular box of fixed area 
A. The expectation value of the observable R = L x /L y was shown to be very different depending 
on whether matter was present or not, and depending on the boundary conditions imposed on 
the matter. In particular, in the presence of a finite density of particles and Dirichlet boundary 
conditions, it was shown that (R) can approach unity and the fluctuations (AR) 2 can vanish 
in the large area limit. This scenario can thus be argued to give rise to a classical, stable, and 
isotropic geometry. A similar effect can arise also in the case of quantum fields if the divergent 
but regularized Casimir energy is included in the analysis. 

The purpose of the calculations presented in this paper is to exemplify that and illustrate 
how matter can affect a dynamical geometry. One technical assumption that was used regarded 



the form of the spectrum for the lengths L x and L y . Although the spectrum was chosen only for 
convenience and not according to a deep physical principle, it should be possible to formulate 
precisely the conditions that must be satisfied by the physical length spectrum for the results 
presented to continue to hold. This kind of analysis, however, goes beyond the scope of this 
paper. Nonetheless, it is worth pointing out that the spectrum used is a natural one from the 
point of view of "atomistic" models of spacetime models in which large-scale geometry is thought 
of being composed of a large number of elementary building blocks. 

A natural application of the presented ideas is in the context of "atomistic" models of 
spacetime. In this context, conventional approaches to studying matter on dynamical geometries 
(see e.g. [2]) have limited use because the notion of geometry is only emergent, but the general 
mechanism described in this paper can be applied. Consider for example the class of models 
called graphity [3,4]: the basic version of these models does not generate extended manifold- 
like graphs, but modified versions that include a homogeneity requirement do [3]. It has been 
suggested that matter degrees of freedom might provide such a homogeneity requirement via a 
mechanism similar to the one described here. Similarly, it is possible that matter might play an 
important role in other models of emergent spacetime such as causal sets [10, 11] or group field 
theory [12]. 

Another application of the ideas presented here is in cosmology: instead of putting 
homogeneity and isotropy as inputs in an ansatz for a cosmological spacetime solution, it may 
be interesting to see them arising from the statistical mechanics of matter fields. 
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